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Abstract — This paper considers feedback load reduction for 
multiuser multiple input multiple output (MIMO) broadcast 
channel where the users' channel distributions are not homo- 
geneous. A cluster-based feedback scheme is proposed such that 
the range of possible signal-to-noise ratio (SNR) of the users are 
divided into several clusters according to the order statistics of 
the users' SNRs. Each cluster has a corresponding threshold, 
and the users compare their measured instantaneous SNRs with 
the thresholds to determine whether and how many bits they 
should use to feed back their instantaneous SNRs. If a user's 
instantaneous SNR is lower than a certain threshold, the user 
does not feed back. Feedback load reduction is thus achieved. 
For a given number of clusters, the sum rate loss using the 
cluster-based feedback scheme is investigated. Then the minimum 
number of clusters given a maximum tolerable sum rate loss is 
derived. Through simulations, it is shown that, when the number 
of users is large, full multiuser diversity can be achieved by 
the proposed feedback scheme, which is more efficient than the 
conventional schemes. 

I. Introduction 

Multiuser diversity can significantly improve system 
throughput when the users suffer channel fluctuations. The 
performance gain of multiuser diversity grows with the number 
of users when the scheduler performs maximum throughput 
scheduling |JJ. For the broadcast channel, the duty paper 
coding (DPC) Q achieves a sum rate which was shown in 
O to have a growth rate M \og\og{K), where M is the 
number of transmit antennas at the base station (BS) and K 
is the number of users in the system. However, this scheme 
has high complexity in encoding/decoding and is difficult to 
be implemented. Therefore, a suboptimal and low-complexity 
zero-forcing (ZF) beamforming technique was proposed in 
ill] which also achieves the optimal growth rate of the sum 
rate. The results both in DPC and ZF schemes were based 
on the full channel state information (CSI) assumption at the 
transmitter and thus the users are required to feedback perfect 
CSI to the BS. Although the optimal sum rate can be achieved, 
the feedback load will increase linearly with the number of 
users. 

Various approaches have been proposed to limit the amount 
of feedback load and investigate sum rate loss. In 0161, the 
quantized channel direction information (CDI) was used to 
characterize the sum rate loss when the ZF technique is con- 
sidered. It was shown that the number of feedback bits of each 
user needs to be increased linearly with the transmission power 



to maintain a constant sum rate loss. Another low-feedback- 
rate and practical scheme, orthogonal random beamforming 
(ORB) was proposed in Q- In ORB, each user only feeds 
back the CSI and the beam index of its favorite beam to the 
BS. Therefore, the total amount of feedback can be reduced 
from MK CSI values to K CSI values. Besides, the sum rate 
loss is negligible when the number of users is large |8|. In 
an effort to further reduce feedback load, a threshold based 
mechanism was proposed in ||9l such that a user does not feed 
back when its CSI is below the threshold. In that work, the 
design of threshold does not take the scheduling algorithm 
into account. In lITOl . multiple thresholds were proposed. The 
design of multiple thresholds was based on the order statistics 
of the signal-to-interference-noise-ratio (SINR) assuming that 
greedy sum rate scheduling was performed in a homogeneous 
network where the users' channel gain distributions are the 
same. Exhaustive search was used to obtain the thresholds, 
which resulted in high computational complexity. 

In this work, the closed form solution of the multiple 
thresholds in [lOj is found, and a more realistic assumption 
on the channel distributions of the users is considered. We 
assume that the distributions of the users' signal-to-noise 
rations (SNRs) are independent and non-identical. Since the 
computational complexity is too high to consider all the users' 
non-identical channel statistics, the statistics of the users' 
channels are divided into multiple clusters. The statistics in 
each cluster are used to calculate the corresponding threshold. 
Each cluster corresponds to an SNR range and is quantized 
with a few bits for differentiating the users' SNR levels falling 
in the same cluster In addition, the sum rate loss due to setting 
thresholds is investigated. The performance of the proposed 
scheme is compared with that of the conventional feedback 
scheme and single threshold feedback scheme [9] in terms of 
sum rate, feedback load and efficiency. Through simulations, 
it is shown that the proposed cluster-based feedback scheme 
is more efficient than conventional feedback schemes and 
achieves higher sum rate than the single threshold feedback 
scheme. 

II. System Model 

The multiuser multiple-input multiple-out (MIMO) down- 
link system is considered. The BS is equipped with AI anten- 
nas and there are K users, each having N receive antennas. 



According to the ORB strategy for multiuser transmission, 
the BS uses a precoding matrix W = [wi, W2, . . . , wm] 
to simuhaneously transmit signals, where £ C*'^^^,^ = 
1,2, . . . , M, are random orthogonal vectors generated from 
isotropic distribution [11]. Let s = J2m '^mSm be the M y.1 
vector of the transmitted signal, where s„i is the mth trans- 
mitted symbol. The received signal for user k is given by 



Yfc = Hj..Ws + rife 



(1) 



where denotes the channel matrix between the BS and 
user k. The elements of the channel matrix are assumed 
to be independent identically distributed (i.i.d.) complex Gaus- 
sian random variables with zero mean and variance aj,. The 
noise term for user k, denoted rifc, is modeled as i.i.d. zero 
mean complex Gaussian with covariance matrix E[nfen^] = 
<t|^I, Vfc, where E denotes the expectation operation and (■)^ 
represents the transpose conjugate. 

User k uses the ZF receiver [il2i to perform channel 
inversion to the received signal y^. Thus, the received signal 
after ZF receiver is given by 



(HfcW)Vfc =s + (HfcW)^nfc 



(2) 



where = (H^H)"iH^ is the pseudo -inverse of H. Under 
the equal power assumption, let the transmit energy of each 
antenna be jj, where P is the total transmit power at the BS. 
The SNR for the /cth user at the mth spatial channel is denoted 

by Xm,k 



P/M 



1,2,...,M(3) 



where [A]„j denotes the mth diagonal element of matrix 
A. Assuming N > M, it is well known that Xm.k is a 
chi-square random variable with 2{N — M + 1) degrees of 
freedom ||T3lllT41 . Then the probability density function (PDF) 
of Xm^fcjVm, can be expressed as 



fk{x) 



(4) 



where \k = ^p°^" . For simplicity, we drop the index m, denote 
Xm,k as Xk, and restrict our analysis for the case N = M. 
Therefore, the distribution of Xk becomes the exponential 
distribution with parameter Xk- 

Let , X(2) , • • ■ I -^{K)^ be the order statistics of indepen- 
dent continuous random variables Xi,X2, ■ ■ ■ , Xx, with the 
PDF (HI in decreasing order, i.e., > X(2) > • • • > -^{k)- 
The cumulative distribution function (CDF) of the largest order 
statistics can be shown as 



K 



XiX\ 



and its corresponding PDF can be expressed as 



(5) 



(6) 




Fig. 1. Cluster-based Model. 



where J^t denotes the summation over all K\ permutations 
{ii, 12, ■ ■ ■ , Ik) of (1,2,..., K). Applying the maximum sum 
rate scheduling algorithm, the sum rate of the system can be 
written as follows: 



R = 




max Xra.k 
l<k<K 



\0g2{l+x)fxnAx)dx. 



(7) 



In order to achieve the sum rate described in (|7}, each user 
should feed back the index of the precoding vector on which 
it sees the highest SNR, and the corresponding SNR value. 

III. Cluster-based Feedback Model 

The transmission and feedback procedure can be described 
as follows. From the previous channel condition or location 
information feedbacks from the users, the BS derives the users' 
mean SNR. The BS then groups the users' mean SNRs into 
Nc clusters according to their magnitudes. The mean SNRs in 
each cluster are similar in quantity and are used to derive one 
SNR threshold, denoted rc,i for cluster i. Note that derivation 
of the users' mean SNRs and the cluster thresholds is done 
periodically. Derivation of the users' mean SNRs does not 
need to be very accurate. The BS broadcasts periodically the 
threshold set {vc^i > rc,2 > ■ • • > ^c.tVc} to the users. 

At every feedback instant, each user compares its measured 
instantaneous SNR with the thresholds to determine which 
cluster its instantaneous SNR belongs to, and feeds back to 
the BS the cluster index and the quantization bits of the 
instantaneous SNR for this cluster When a user's SNR is 
smaller than rc.N^, the user does not feedback to the BS. The 
proposed procedure takes a little downlink bandwidth for the 
BS to periodically broadcast the threshold set, in exchange of 
reduction of the uplink feedback bandwidth. 

To derive the thresholds, the overall statistics of users 
are divided into multiple clusters according to the mean 
SNRs T- , T- , • • • , -T- ■ The means of the random variables 
Xi , X2 , ■ ■ ■ , Xk are ranked with decreasing order as > 
-p!— > ■ • • > — !— and uniformly divided into Nc clus- 

A(2) ~ ~ '^(K) 

ters with size L = K/Nr. as \ -r^, • • • , • • , 

' L A(i) A(2) A(i) J 



< -r , , • • • , f ■ To simplify the nota- 

l.^(i(Nc-l) + l) -*'(I.(iVc-l) + l) -^(K) J 

tion, we let random variable represent the mth statistic 
in the nth cluster Then, Y"^ is exponentially distributed with 



mean value 



The random variables in cluster n 



are {F^", Kf , • • • , F/'}, and let Y^l^ > Y^^^ > ■ ■ ■ > Y^l^ 
be the order statistics of them. For a measured instantaneous 
SNR Hi at user i, we say that the rank of i/i in cluster n 
is d if Y^l > >Y^ly> y, > . . . > Y(i,y In 
" V ' " V ' 

{d — 1) variables (L — d) variables 

the following, we will discuss how to derive the threshold 
in each cluster for heterogenous and homogeneous channel 
distributions to reduce the feedback load. 

A. Heterogenous Case 

1 ) Cluster-based Type-I: With maximum sum rate schedul- 
ing, a user will be scheduled when its SNR on a particular 
precoding vector is the highest among the users. On the other 
hand, if a user has a low SNR, it is unlikely to be scheduled. 
For this user, feeding back CSI is wasteful of the uplink radio 
resource. The threshold of each cluster is designed according 
to the probability of a user's measured instantaneous SNR 
being a particular rank in that cluster Let „ (r) denote the 
probability that user m is ranked n among the L users in 
cluster i when its instantaneous SNR is r. 

= >Y^>--> Yl_^ > r > Yl^^ • • • > 



(n — 1) variables 



(L — n) variables 



b=n+l 



where denotes the summation over all permuta- 

tions (tl, ■ • • ,t„_l,t„+l, ■■■ ,tL) of (1,2, ■ • • ,TO - 1,TO + 

1, • • • , i) in the cluster i. For example, when the instantaneous 
SNR of user t in cluster 1 is infinity, the rank of user t 
among the L users in cluster 1 is one with probability one, 
i.e., Pt^i(oo) = P{Yt^ = = oo} = 1. Being valid 

conditional probabilities, the ni^Y^ satisfy 



£p4,„(r)-l, m-l,2,..-,L. 



(8) 



The most probable rank of user m in cluster i when its 
instantaneous SNR is r can be obtained by 



P* (r). 

712,71 \ J 



Let Qi 



arg max ^ 

ne{l,2,...,L} 



be defined such that PI 



(9) 



{Qm.ri) 



P'mn+i {Qnni}- ^^'^ ^ic Seen that when the instantaneous 
SNR of user m falls in the range [QJ^ „, (35„ the most 

probable rank of user m in cluster i is n. For maximum sum 
rate scheduling, we let only the users who are most likely to 



be rank one in each cluster to feedback. Therefore, a threshold 
for each cluster i is set as 

r,,, = max{Ql_i, Q^^, • • • , Ql^,} i = 1, 2, • • • , iV, (10) 

which satisfies the following inequality 

rc,i > rc,2 > ■ ■ ■ > rc,N,- (H) 

To reduce the computational complexity, another simple and 
low complexity method is proposed in the next section. 

2) Cluster-based Type-II: We approximate by assuming 
that the random variables Y^,m = 1,2,- •• ,L in cluster i 
have the same exponential distribution with mean obtained 
by 



fJ-c 



i~l)+m) 



(12) 



Then, P,*^ „(r) is the same for all users in a cluster, and can 
be obtained by 

p, / N (-^-l)!exp(-/Xe,r)"~^(l-exp(-^e.O)^~" 
"^"^ ' (n-l)!(L-n)! ■ ^ ' 

With this approximation, the closed-form solution of Tc^i can 
be derived as 



1 



\n.L,i = 1,2, 



(14) 



B. Homogeneous Case 

The homogeneous case can be viewed as a special case 
of the cluster-based type-II scenario using only one cluster 
{Nc = l,K = L). The mean values of the random variables 
Y^, TO = 1, 2, • • • , if are the same, i.e., l/A = l/Ai = -- - = 
1/Xk- As in ifTOl , multiple thresholds are set according to the 
most probable rank as 



(15) 



For all cases, at each feedback instant, when a user's 
instantaneous SNR is greater than the smallest threshold Tc^at^, 
the user needs to feed back to the BS using Be ~ log2 [A^d 
bits to indicate which region (between two adjacent thresholds) 
its instantaneous SNR belongs to. In addition, in order to 
differentiate the users SNR in the same region i, the region i 
which is represented by £j is further quantized with bi bits. 
Therefore, the feedback bits of each user include two parts: 
one is the region index with Bq bits and the other is the 
quantization bits for that region. 

IV. Analysis of Sum Rate Loss 

A natural question to ask is how many clusters (thresholds) 
should be set. Obviously, if many clusters are used, the number 
of region index bits Be is increased. On the other hand, when 
a small number of clusters are used, the smallest threshold 
Vc^Na becomes large. Then, the sum rate loss caused by 
the threshold rcjv^ may increase. Therefore, an appropriate 
number of clusters is important for the design. We now 



analyze the sum rate loss of the system caused by the smallest 
threshold Tcjv ■ We define a random variable Zi as follows: 



0, 



(16) 



According to the maximum sum rate criterion, the exact SNR 
loss for the system is = maxjZi, Z2, • • • ,Zk}- The 

probability of the rate loss event can be obtained by 



Pl = G {0,ra^Nj} = [](! -cxp(-A,r,,wJ). (17) 



i=l 



Therefore, the sum rate loss can be expressed as 

Ai?- A/S{log2(l + Z(i))}Pi. 



V. Bit Allocation and Feedback Load Analysis 

A. Bit Allocation 

Let Tc.o = 00. Assume that the SNR region i [rc,i,rc,i-i), 
denoted = 1,2,- ■■ ,Nc, is quantized with bi bits. In 
region £i, the quantization levels using bi bits are expressed by 
ql,t = 1,2, ■■■ , 2^\ obtained by a pdf quantizer Thus, 
each level will occur with the same probability. Under the 
per user average feedback load constraint Ck for user k, the 
available bits will be assigned to the regions to maximize sum 
rate. Let Pc, = P{X(i) G = 1,2,--- ,iVc- The bit 

allocation problem can be described as follows: 

; iog(i + q\) 



(18) 



Theorem 1. 4751/ Let the means and variances of the random 
variables Zi, Z2, ■ ■ ■ , Zk be ^ — {^i, ^,2, - ■ ■ , IJ'k) ond cr^ = 
(ctj , erf , • ■ • , (^k), respectively. The closed form upper bounds 
on the expected value of the largest order statistic is: 

K 



E{Z^,)} < Y.{ 



IM + \J (Mi ~TY + erf 



where T = maxi<j<K{//j + j-^==aj}. 

Applying Jensen's inequality and the upper bound of 
E{Z^i-^}, the sum rate loss on a certain beam can be bounded 
by 

A 7? 

— = E{\0g2{l + Z^,))}PL<\0g2{l+E{Z^^,)})PL 

< log^il + fiz^^^)PL. (20) 

Using ( |20] i. for a given tolerable sum rate loss ARu, the 
minimum number of clusters (thresholds) required can be 
determined. 

Proposition 1. When the total number of users approaches 
infinity, the sum rate loss caused by the smallest finite thresh- 
old rc^Na negligible. Thus, the full multiuser diversity can 
be achieved. 

Proof When the total number of users K approaches to 
infinity, lim P^ = 0. In addition, Z(i) is bounded by Vc n^- 

K—>-oo ' " 

The sum rate loss on a certain beam becomes zero 
Ai? 

lim — = lim i?{log2(l + Z(i))}Fi 

X— i-oo Ivl 

< lim log2(l + S{Z(i)})FL 



< lim log2(l + rc,jvJ-PL 



(21) 



max 

(hi,h2,--- ,bNa) 



s.t 



E 



Pu.M <Ck,k = 1,2,- ■■ ,K, (22) 



where P^^ ^ = P{Xk G Cj} = e 



(2 - K)T 



(19) 



g-Afcre.,_i^ The 
closed-form integer solution for (122) does not exist. We use 
exhaustive search to find the optimal bit allocation set to 
maximize the sum rate. 

B. Feedback Load Analysis 

When a user's instantaneous SNR is smaller than the 
smallest threshold rc^N^, the user does not need to feedback. 
Thus, the feedback probability for user k is e^^'''"'^ "<=. The 
average total feedback load can be expressed as 



VI. Simulation Results 



(23) 



In this section, we show the sum rate and feedback load 
performance for the proposed feedback scheme. The BS is 
equipped with M — A antennas and the total number of users 
is i^T = 10 ^ 100. In order to perform ZF beamforming, we 
let the number of receive antennas = 4. The total transmit 
power P is lOW, while the additive white Gaussian noise 
power at the receivers a'j^ is IW. Note that these numbers 
are selected only for illustration purpose. The elements of 
the channel matrix Hi for the ith user are assumed to be 
i.i.d. complex Gaussian distribution with zero and variance af, 
where af are drawn uniformly from the interval [0,1] to model 
heterogeneous channels. The number of clusters in type-I and 
type-II schemes is set to four (thus Be = 2) according to the 
tolerable sum rate loss /S.Rjj = 10^^ bps/Hz. The feedback 
load constraint of user fc is = 0.8, fc = 1, 2, • • • ,K.\n the 
simulation, the proposed feedback schemes are compared with 
the conventional scheme and the single-threshold schemes 
proposed in |9|. In the conventional feedback scheme, no 
matter what the instantaneous SNR is, it is always quantized 
with 3 bits. In the single-threshold scheme, the region {rth, 00) 
of the SNR is quantized using 3 bits. The single threshold rth 
is established according to the scheduling outage probability 

Pout- 




Number of users 



Fig. 2. Sum rate comparison between different feedback schemes. 
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Fig. 4. The efficiencies of different feedback schemes. 
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Fig. 3. Feedback load comparison between different feedback schemes. 

In Fig. 121 the type-I feedback scheme achieves the highest 
sum rate, and the low-complexity type-II scheme achieves 
almost the same rate as the type-I scheme. As sho'wn in Fig. |3] 
the feedback load of the conventional scheme increases lin- 
early with the number of users. Using our proposed schemes, 
the total number of feedback bits can be dramatically reduced. 
Overall, the proposed schemes use fewer bits to achieve higher 
sum rate than the conventional scheme. The single-threshold 
scheme has lower feedback load, but suffers significantly in 
the sum rate performance. 

In Fig. m we plot the sum rate vs. the total feedback load 
as an indication of the efficiency. The type-I scheme is the 
most efficient (i.e., making best use of the feedback bits), but 
has high computational complexity. The low-complexity type- 
II scheme not only achieves high sum rate but also reduces 
the feedback load significantly. The single-threshold feedback 
scheme only achieves a sum rate of about ll{bps/ Hz) with 
a small feedback load. 

VII. Conclusion 

In this paper, we investigated the feedback load reduction 
problem in multiuser MIMO broadcast system. We proposed a 



cluster-based feedback scheme to reduce the feedback load in 
heterogenous and homogeneous channels. The bit allocation 
problem for the multiple-cluster feedback scheme was also 
discussed. The simulation results showed that, compared to the 
existing feedback schemes, the cluster-based feedback scheme 
can make the best use of the feedback bits to achieve good 
feedback load reduction while maintaining good sum rate 
performance. 
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